Approximate analytical description of the nonaffine response of amorphous solids 
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An approximation scheme for model disordered solids is proposed which leads to the fully analyt- 
ical evaluation of the elastic constants under explicit account of the inhomogeneity (nonaffinity) of 
the atomic displacements. The theory is in quantitative agreement with simulations for central-force 
systems and predicts the vanishing of the shear modulus at the isostatic point with the linear law 
fl ~ [z — 2d), where z is the coordination number. The vanishing of rigidity at the isostatic point is 
shown to be the consequence of the canceling out of positive affine and negative nonaffine terms. 

PACS numbers: 46.25.-y, 64.60.aq, 45.70.-n 



I. INTRODUCTION 

Disordered or amorphous solids represent a great part 
of ordinary matter (e.g., glass), including biological mat- 
ter (e.g., cytoskeletal networks) [THS]- Yet, the relation- 
ship between rigidity and disorder has remained elusive 
and no theory has hitherto proved able to correctly de- 
scribe their elastic constants. The rigidity of disordered 
solids is also intimately related to the fundamentally un- 
solved problem of the glass transition |4J. Elastic rigid- 
ity in supercooled liquids emerges from the fluid state 
at the glass transition without any detectable lowering 
of the symmetry (apart from translational and replica 
symmetry-breaking), as opposed to what happens in "or- 
dinary" liquid-solid phase transitions 5] . In recent years 
it has been recognized that the intrinsic " softness" of dis- 
ordered solids is related to the nonaffinity of the atomic 
displacements [TJ HI |6j [7] : the atoms in a strained disor- 
dered solid are not displaced proportionally to the global 
strain. The calculation of the contribution to rigidity due 
to such nonaffine displacements poses formidable difficul- 
ties because it requires the analytical knowledge of the 
cigenmodes of the dynamical or Hessian matrix of the 
system [5], which is a sparse random matrix. The 
simplest disordered solids where nonaffinity is supposed 
to play a significant role are those with nearest-neighbor 
central- force interactions [HI Ell- In these systems, it is 
well-known that the shear modulus \x vanishes at the iso- 
static point where each particle has an average number of 
nearest neighbors (mechanical contacts) z — 2d [TO] . 
Although it seems reasonable from constraint-counting 
arguments that rigidity is lost when z = 2d [TT], the 
linear law ji ~ (z — 2d) observed in simulation studies 
of completely disordered solids [TO] has remained un- 
explained and the physics behind it represents a long- 
standing problem [TO] where the role of nonaffinity is 
yet unclear. Further, the well-documented inadequacy 
of affine theories to describe the elasticity and transport 
properties of amorphous materials calls for an improved 
theory beyond the affine approximation |13) . In this Let- 
ter, we propose an approximation scheme which gives a 
well-defined deterministic limit for the nonaffine contri- 



butions to the elastic constants. This leads to a fully 
analytical description of the elastic constants which ac- 
counts for the microscopic nonaffinity of the atomic dis- 
placements. 



II. FORMALISM 

In the following, Roman indices are used to label atoms 
while Greek indices are used to label Cartesian compo- 
nents. The summation convention over repeated indices 
holds throughout for Greek indices. Bold characters de- 
note vectors in diV-dimensional space (7V=total number 
of atoms). We closely follow the notation of Lemaitre 
and Maloney |5] and we start our analysis from the def- 
inition of a Bravais cell for the disordered lattice. The 
cell is described by three Bravais vectors and thus by a 
matrix h. The potential depends on the particle position 
Tj and on the shape of the cell which enforces bound- 
ary conditions. Macroscopically imposed deformations 
of the cell are described by changes in the Bravais vec- 
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tors through a linear map F — h h and the relation 

h = Fh, We denote quantities in the reference frame, as 
well as quantities that are measured with respect to the 
reference frame, with a circle. Accordingly, the unit cell 
in the reference configuration is given through the ma- 
trix h. In the language of continua F is the deformation 
gradient tensor. The position of atom i after an affine 
deformation is given by 

r? = F aP RP (1) 

As a unique exception to the ring notation, we denote 
with J?j the position of the atoms in the reference cell. 
This relation illustrates the definition of affine displace- 
ment as atom i is displaced proportionally to the external 
deformation. The position of the atom after a strain is 
denoted by r" (F) and differs from Eq. |T|) if the nonaffine 
displacement is not zero. It is useful to introduce the 
particle position for an atom which undergoes both 
affine and nonaffine displacements, defined by 



r?(F) = F^r^F). 



(2) 



2 



Remark also that R t = £,(0). Hence, for an affine dis- 
placement the position vector is kept fixed while the 
new position is determined by F_. Any additional (non- 
affine) displacement is thus parameterized for a given 
strain in terms of r$(£). With these definitions we can 
express the potential U ({z^}, £) in the coordinates of the 
reference frame as W({r i },.F) which is defined by 



(3) 



The elastic constants are defined by 
1 d 2 U 



(6) 
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In order to account for the nonaffine relaxation in the cal- 
culation of the elastic constants, the derivatives in Eq.Q 
have to be taken along a trajectory of (locally) minimum 
energy. Following Lemaitre and Maloney [8], we obtain 



It is convenient to introduce the Cauchy-Green strain 
tensor 77 = ^(F 1 F_ — F) to describe the deformations 

since the elastic constants are defined in terms of second 
derivatives in r\. The deformation is completely described 
by this tensor since the total internal energy of the solid 
being deformed can be expressed as W({|zVj|}) ; i- e - a s a 
functional of the set {r^-} of relative distances between 
atoms Tij = \r_ij \ = \r_i ~ Tj|> and 77 describes affine trans- 
formations in terms of relative interatomic distances ac- 
cording to \ERij? = lEij? + />',,'/ />',,■ 

A homogeneous strain of the cell F will first bring each 

atom to its affine position rf = F a pR^. In this affine 
position the total force acting on atom i is in general 
not zero since the neighboring atoms may exert a non- 
vanishing force-field due to their affine motion. This is 
especially true for disordered solids where the nearest 
neighbors are placed at random around atom i so that 
they transmit unbalanced forces to i (in an ordered lat- 
tice the transmitted forces balance by symmetry such 
that this effect is often negligible). It is in response to 
these virtual forces that the atom undergoes an addi- 
tional motion after it has been displaced affinely such 
that the energy released in the process reestablishes (lo- 
cal) mechanical equilibrium. Hence the system under 
reversible strain evolves adiabatically along a trajectory 
r(r)) that minimizes the mechanical energy for a given 

strain r). If we denote by v ^ the derivative with respect 

to adiabatic changes of the strain under the constraint of 
mechanical equilibrium, one obtains an equation of mo- 
tion of the nonaffine displacement by differentiating the 
force f" = — evaluated in the true position (where it 
vanishes). In the limit 77 — >■ the equation reads 



(4) 
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where we have introduced the Hessian Hff and the affine 
force field S° given respectively by 



H 



a& = d 2 U 
« dffdrf 
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(7) 

where it is evident that the true elastic constant is given 
by the affine (Born-Huang) elastic constant C4 cor- 
rected by the nonaffine term —C^^ % . Following Lemaitre 
and Maloney [8] , and using Eq. (4)inEq. (7), one derives 
the following expression for the nonaffine correction, 



s-iNA 



(8) 



The last inequality in Eq. (8) is justified in view of the 
Hessian matrix being semi-positive definite at mechanical 
equilibrium. Hence it follows that the correction due to 
the nonaffine relaxation 
negative contribution to the total rigidity [H] 



< 0, necessarily gives a 



III. APPROXIMATION SCHEME 

A. The Cauchy bonded-network model 

Let us consider the disordered Cauchy solid, defined 
by the following properties pQ: (i) atoms interact pair- 
wise and only with their nearest neighbors; (ii) the in- 
teraction potential is a central-force harmonic potential; 
(Hi) the reference state is unstressed, i.e all springs (in- 
teratomic bonds) are relaxed in the minimum of the 
harmonic well; (iv) disorder is spatially decorrelated. 
The equivalence with a random network of harmonic 
springs is evident [H [14] . Hence, the total free energy is 
given by W({ry}) = J2(ij) v ij( r ij) where the sum runs 
over all pairs of nearest- neighbors (ij). The pair in- 
teraction potential is given by the harmonic potential 
V(rij) — § (rij — Ro) 2 - k is the atomic force constant 
and Ro is the interatomic distance at rest in the refer- 
ence frame. Under these conditions the Hessian matrix 
becomes 

H if = S v MHsnf s n% a - (1 - S tj )Kc ij nf J n^ (9) 
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where we used the identity d/d^ 



jd/drij, with 



TLij = Lij/ r ij- Further, Cjj is the (random) occupancy 
matrix with Cjj = 1 if i and j are nearest neighbors and 
Cij = otherwise. Cjj is a matrix where each row and 
each column have on average z elements equal to 1 dis- 
tributed randomly under the constraint that the matrix 
be symmetric. Using this form of the Hessian one obtains 
the affine part of the elastic constant as 



/~iA RqK V A L £ K X 

ij 



(10) 



which is the well-known Born-Huang formula [TJ UJ 
|S]. Further, we also obtain a microscopic expres- 
sion for the affine force field from Eq.|5l) as E? 



J2i R ij> 



, . . , t . .CijUfjU^nfj . We can now turn to the non- 
affine part of the elastic stiffness, Cfl„. ■ The Hessian is 
a dN x dN symmetric semi-positive definite matrix with 
d eigenvalues equal to zero which are due to the global 
translational invariance of the solid. Eq. ^ can be solved 
by normal mode decomposition which leads to [5] 
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(11) 



where y fc are the eigenvectors of the Hessian (which are 
orthogonal since the Hessian is symmetric), A^ the corre- 
sponding eigenvalues and (, ) denotes the normal scalar 
product on M. dN . In the next section, we shall evalu- 
ate the deterministic limit of Eq. (11), which is a self- 
averaging quantity [S]. 



B. The approximate Hessian and the affine field 
projection on its eigenmodes 



implies sacrificing some details of the vibrational spec- 
trum. This problem is addressed in section IV. B and 
IV. C where we show what details are lost and we study 
the validity and limitations of the approximation. 



In d = 3 it is — (cos (j>ij 
and the pair of angles 
orientation of the bond (ij) 



sin O^j , cos @>ij 



ij and 9ij univocally specifies the 



factors in the Hessian, nf„nf. 



The orientation-dependent 
in Eq. (9), for a large 
system with uncorrelated isotropic disorder (where every 
bond can take any orientation in the solid angle with the 
same probability 1 /in) , can be replaced with its isotropic 
(angular) average, i.e. n^n^- 5 a p/d. Within this ap- 
proximation, the Hessian becomes 



(i 



Sat 



(12) 



Remark that this is still a sparse random matrix because 
of the positional disorder in the random coefficients Cij. 
According to Eq.(12|, let us define H = H_® I where 
I is the d x d identity matrix (which represents 6 a p) 

and H is the matrix which multiplies S a p in Eq.(12). 
Denoting with {a q } q= \..N the set of eigenvectors of H_, 
which is an orthonormal basis (ONB) of M. N , and with 
{&i}i=i..d the standard Cartesian basis of K d , it follows 
that (H_ ® I)(a<E) e) = A (a ® e) and thus the dN dimen- 
sional set {a q e l } q =i..Nj=i..d is an ONB of eigenvectors 



of as given by Eq.(12|. This allows us to write (with 
v = ae ; for some a <E \Q q }q=i..N)'- 



A rigorous derivation requires one to first determine 
the eigenmodes v fc of the Hessian matrix given by Eq.(9) 
in order to calculate the projection on them of the affine 
fields in Eq.(ll). Thereafter the average over the dis- 
order is taken to get the thermodynamic limit of C^^.. 
However, the Hessian is a random matrix and both y fe 
and Afc depend on the realization of disorder. Also, be- 
ing the Hessian sparse, there are no analytical forms even 
for its statistical spectral distributions. Nevertheless, the 
deterministic limit of Eq.(ll) can be calculated analyti- 



cally within the following approximation that we propose 
here. 

Our approximation consists in performing a disorder- 
average of the orientation-dependent part of the Hessian 
first and then use the result to calculate the eigenmodes, 
their inner products with the affine fields and finally the 
nonaffine correction. Inverting the sequence of " calculat- 
ing" and "averaging" is sometimes referred to as an ef- 
fective medium approximation and is not at all unusual 
in dealing with disordered systems [T5] since it is often 
the only strategy to keep the treatment analytical. Us- 
ing some averaged form of the Hessian matrix necessarily 




r sr' s' 



(E^y) (H kx ,y) - ( ^a r %^ ; \ «>-,,,,, > < 

\^{(X r (X r ' C rs C r ' s ' ) 

x (n rs n rs n^ s n r / S /Ti r / s /7i r , s ,)} 

With our isotropic approximation, we replace the 
orientation-dependent terms with their isotropic angular- 
averaged values which gives n\. g n^ g nf g n r , s ,n^, s ,n^, s , = 
(6 rr >5 ss , - S rs 'S sr ')-B L ^ KX where the B L ^ KX are geomet- 
ric coefficients resulting from the angular average. For 
d = 3 and d = 2 they are as follows 
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Substituting in Eq.(ll| we obtain 
X I ^ ^ Gy c rs c r _ 



/ v Oj r Cl s C rs C sr I 
r s / 



— K 2 Rl Bi^ kx — ^ a r a s H r 



(15) 



where we used that c 2 rs = c rs c sr = c rs and the identities 



X]rs ^r^s^-rs — X]rs a r a s [(^^/j Crj*)&rs 

c rs (l - S rs )] = ^J2^s a r a sHrs- Recalling that 
J2^H rs a s = Xa r , we obtain (S^,v fc ) (H KX ,v fc ) = 



IV. RESULTS AND DISCUSSION 
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FIG. 1: (color online). The theoretical predictions for the 
shear modulus, Eq. (17), and the simulation results of 
Ref. 1 12], No fitting parameter is used in the comparison, 
re = 1, Rq = 1. 



Elastic moduli 



Hence, we have shown that within the isotropic ap- 



proximation of the Hessian, Eq.(12|, the nonaffine part 



of the elastic stiffness, Eq.(ll I, has the following thermo- 
dynamic limit 



N d 



V , , , A q 
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(16) 



The affine part of the elastic constants for the disordered 
Cauchy solid can b e o btained by performing the disor- 
der average of Eq.(10l, (C^ ), where (.) denotes the 
angular average, and we always use the isotropic distri- 
bution of the bond orientations. In d = 3 we thus obtain 



/i = {C^y ) = T^yKzRQ, for the shear modulus, and 

K A = l((C* xx ) + 2(C4 W )) = ^kzRI for the bulk 
modulus. Therefore, using these affine moduli together 
with the coefficients of Eq.(14) and with Eq.(16) we de- 



rive expressions for the shear and bulk modulus of the 
d = 3 disordered Cauchy solid, respectively as 



K = K / 



For d = 2 we obtain 
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K = K A - K NA 
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1 N 
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Generalizing this result to arbitrary space dimensions 
gives the following scaling for the moduli in d dimensions 



PL ~ K ~ (z - 2d) 



(19) 



The predictions of Eq.(17) for the shear modulus, with- 



out fitting parameters, can be compared with the simula- 
tions of Ref. [H] of d = 3 disordered packings of (monodis- 
perse) compressible spheres interacting via harmonic re- 
pulsion in Fig. ([I]). In the simulations the spheres, at 
T = 0, are slightly compressed to packing fractions 4> just 
above the so called jamming point at (f>j — 0.64 which is 
also an isostatic point with zj = 2d [HI H2] • Further, the 
jamming point is a zero-stress point |12j . and the effect 
of stress on the global rigidity is therefore small. How- 
ever, it seems from this comparison, and from our results, 
that stresses, in general, are not likely to affect the qual- 
itative behavior of the global rigidity as they play no role 
in arriving at the fundamental scaling law Eq.(18). 

Finally, we note that the vanishing of K at the isostatic 
point predicted by our theory does not agree with the 
scaling of K observed in soft-sphere packings where it 
remains finite at zj = 2d [T^l E]. It agrees however 
with the behavior of random networks where K vanishes 
linearly at zj = 2d [14] . The reason for this might be 
tentatively identified with the fact that in our theory, just 
like in networks, excluded volume effects are irrelevant, 
whereas they are important in packings |14) . 



B. Vibrational density of states and validity of the 
approximation 

The isotropic Hessian matrix introduced for the calcu- 
lation of the nonaffine contribution to the elastic moduli 
can be used to obtain the density of vibrational states 
(DOS) numerically. Recall that the approximate Hessian 
is given by the following random matrix 



5ij)Ci 3 
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The Hessian is defined by the coefficients Cy which de- 
pend on the realization a and are, therefore, random 
variables. The eigenvalues and thus the eigenvalue dis- 
tribution of a random matrix are also random quantities. 
The explicit calculation of the eigenvalues as functions 
of the matrix elements is not possible. The approach to 
the eigenvalue problem in random matrix theory makes 
use of the self-averaging assumption that the eigenvalue 
distribution becomes deterministic in the limit of an in- 
finite system size [16]. As analytical solutions for the 
eigenvalue distribution in our case are not possible (due 
to the sparseness of the Hessian) , we resort to a numer- 
ical analysis assuming that the self-averaging property 
holds. Therefore, we can define a limiting eigenvalue dis- 
tribution p(X) as follows 



lim (p N (\)) a = P (\) 



N 



Then we have for all realizations a that 



N 



lim p N (X)[a] = p(A) 



Remark that for a finite N the eigenvalues distribution 
is discrete and given by 



p N (X)[a] 



1 N 



where 5 is the Dirac delta function. In the limit N — > 
oo, the set of eigenvalues has infinite elements, and the 
distribution becomes continuous. The factor 1/N is nec- 
essary to normalize the density given that the normal- 
ization condition is J Q p(X) dX = 1. To analyze numer- 
ically the eigenvalue distribution we calculate the eigen- 
values sets {Xi} r for large systems ( N ~ 10000 ). Then 
we create an histogram of the eigenvalues set and we fit it 
with a continuous curve which approximates the limiting 
eigenvalue distribution. The eigenvalue distribution is 
usually described in terms of the vibrational DOS which 
we denote as -D(w), where u is the vibrational frequency. 
The latter is related to the eigenvalue distribution by the 
change of variables 



and its inverse 



Hence, with du = 



1\Jm\ 



dX and dX — 2mu du we get 



that the density distributions of A and ui are related by: 



D(u) = p(muj 2 )2uj 



and 



P (X) = 



D(\ / A ) 



We compare the so obtained DOS from Eq. (12) with 
the DOS from simulations of unstressed harmonic pack- 
ings [17] where <f> — <f> c = 0.1 (<f) c « 0.64 is the jamming 
packing fraction), corresponding to z — 6 ~ 2. The com- 
parison is shown in Fig. 2. While the upper end and the 
main features of the spectrum (width and average) are 
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FIG. 2: (color online). 1: The DOS calculated for a system of 
TV = 10000 particles according to Eq. (1) (solid line) and z = 
8. 2: Simulations for repulsive unstressed harmonic packings 
of [T7] with <f> — cj> c = 0.1 which corresponds to z — 6 ~ 2. 



quite well reproduced, this was somewhat expected be- 
cause sacrificing information about the bond orientations 
does not change the single-particle parameters (of order 
\J n/ui) which dominate the highly (Anderson) localized 
high-w modes. On the other hand, it is now clear from 
this comparison what details of the vibrational spectrum 
went lost in our isotropic approximation of the Hessian, 
Eq. (12): the spectrum calculated with the isotropic 
Hessian is significantly depleted of low-frequency modes 
whose density is significantly underestimated in compar- 
ison with the DOS of the packings. This is probably re- 
lated to the anisotropic character of the local correlations 
between particles in the packings (ultimately related to 
excluded-volume and static effects) which play an impor- 
tant role in the low-frequency modes and that are lost in 
the approximation. This observation hints again at the 
possibility that overall our model describes random net- 
works (where excluded-volume are absent) better than 
sphere packings. The relationship between these obser- 
vations and the excluded-volume effects should be more 
systematically investigated in future studies. This was al- 
ready noted in relation to the bulk modulus prediction of 
our theory which agrees indeed with the random-network 
scaling (K ~ z — 6), though not with the sphere packing 
one, for arguably similar reasons. 

In view of these considerations, it is natural to ask 
why our theory, which seems to underestimate the low-w 
modes, still yields a correct prediction of the shear mod- 
ulus for sphere packings. This question is related to the 
issue of the role played by the low-frequency modes in the 
nonaffine response. While this issue is a very open and 
unsolved one in our current understanding of amorphous 
solids |10j , a tentative, and certainly incomplete, answer 
to this deep question is proposed in the next section. 
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C. On the role of low-frequency modes in the 
nonaffine response 

To assess the relative importance of different regimes 
of the vibrational spectrum in the nonaffine elastic re- 
sponse, it is instructive to rewrite the elastic moduli with 
the nonafhne correction in the continuous frequency do- 
main. According to the nonaffme linear response formal- 
ism of Lemaitre and Maloney Ref . [5] , in the thermody- 
namic limit one has 

«w = «&*> - f D{u Z% x{u) duJ (20) 

Jq miu 

where the correlators on the frequency shells are defined 
by 

r^ KX (u) = ((B tS , Y fc )(S KX ,v fc )) a , fee [ a , ia)+rfu) ] (21) 

The function r t £ KX (w) thus represents the projection of 
the affine fields on the frequency shells and its magni- 
tude gives the importance of the contribution of each fre- 
quency shell to the nonaffine response. From Eq.(20) it is 
evident first of all that in the zero-frequency limit u> — > 
the moduli diverge to minus infinity, i.e. (C,,£ KX ) — > — oo, 
unless either D(u = 0) = or T^ kx (lj) = 0. At the iso- 
static or jamming point of sphere packings one has that 
the DOS develops soft modes with D{lu = 0) 7^ 0. As the 
nonaffine linear formalism is an exact theory, it is then 
strictly necessary that 

Hmr^w) = (22) 

i.e. the zero-frequency modes must not contribute to the 
nonaffine response. This is what has been observed in- 
deed in the numerical simulations of Ref. [8] where, in 
the case of a Lennard- Jones glass, the function T L ^ KX (u>) 
measured in the simulations goes to zero at oj = 0. Fur- 
thermore, in the same simulation study [8], it was found 
that T L £ KX (uj) not only goes to zero at zero frequency, 
but is a monotonically growing function of ui in the en- 
tire domain, such that it has significantly lower values 
at low w than in the middle and upper part of the spec- 
trum where it reaches its maximum value (cfr. Fig. 5 in 
Ref.JS]). Hence, the simulation results of [5j indicate 
that the contribution of low-frequency modes to the non- 
affine response is small whereas the leading contribution 
comes from the high-frequency modes. This observation 
is also in agreement with physical intuition: the source 
of the nonaffine response is given by the projection of the 
affine fields on the eigenmodes which has a higher value 
the more energetic the modes are. 

Based on these observations, one can conclude that 
the low-frequency modes play a relatively minor role in 
the nonaffine response as compared to the high-frequency 
modes. This explains why our theory, which underes- 
timates the low-frequency modes in the case of sphere 
packings, still yields correct predictions for the shear 
modulus in excellent agreement with simulations (Fig.l). 



In the case of the bulk modulus, simulations [5| give 
practically the same behavior for the correlator T L ^ KX (u>) 
as for the shear modulus, with the low-frequency modes 
contributing to the nonaffine response to a minor extent. 
In this case, the failure of the theory in predicting the 
correct scaling for sphere packings (despite being suc- 
cessful for networks) is more likely to be ascribed to the 
geometric attenuation of the random affine fields under 
hydrostatic pressure due to excluded volume, as we spec- 
ulated in section IV. A. However, this hypothesis has to 
be tested in future work by means of ad hoc numerical 
studies as the bulk modulus scaling of packings is a prob- 
lem currently under debate [14 . 



V. CONCLUSION 

We have developed an approximate, fully analytical 
theory of the nonaffine elastic response of amorphous 
solids which explicitly takes into account the nonaffin- 
ity of the atomic displacements. We have applied the 
nonaffine linear formalism in the formulation of Lemaitre 
and Maloney [5] to the so-called Cauchy bonded- network 
model pQ, i.e. to networks of harmonic central- force 
springs. In order to evaluate the nonaffine correction 
to the elastic moduli analytically, an approximation of 
the Hessian matrix has been proposed where the bond 
orientation-dependent factors in the Hessian are replaced 
with their isotropic average (isotropic Hessian). Even 
though the isotropic Hessian has a density of states 
which significantly lacks low-frequency modes in com- 
parison with sphere packings, our approximation yields 
predictions of the shear modulus in excellent quantita- 
tive agreement with simulations of sphere packings |12j . 
The good agreement is explained with the observation, 
supported by simulations in the literature [5], that the 
low-frequency modes, underestimated by our approxi- 
mation, play a relatively minor role in the nonaffine re- 
sponse, which is controlled by the upper part of the vi- 
brational spectrum (that is well reproduced by our the- 
ory). While our approximation is not suited to accurately 
describe transport properties of disordered solids in the 
low-connectivity and low- frequency limits |18j . it seems 
on the other hand successful in accurately describing the 
elastic response to shear of disordered solids. Further- 
more, our theory provides a completely new insight into 
the linear vanishing of shear rigidity at the isostatic point 
(z = 2d) of disordered solids: this happens because the 
nonaffine correction at the isostatic point becomes equal 
in absolute value, but with opposite sign, to the affine 
part of the shear modulus. 
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